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O . ABSTRACT 



o . 

. The H I mass function (HiMF) for galaxies in the local universe is constructed 

G^ 
^. 

p ^. in the 21cm line. The survey consists of two strips covering in total ~ 65 

square degrees of sky, with a depth of cz = 7400 kms~^ and was optimized 



to detect column densities of neutral gas A'hi > lO^^cm ^ {5a). The survey 



O 

! yielded 66 significant extragalactic signals of which approximately 50% are 

cataloged galaxies. No free fioating H I clouds without stars are found. 
^ ! VLA follow-up observations of all signals have been used to obtain better 

measurements of the positions and fiuxes and allow an alternate determination 
of the achieved survey sensitivity. The resulting HiMF has a shallow faint 
end slope {a ~ 1.2), and is consistent with earlier estimates computed for 
the population of optically selected gas rich galaxies. This implies that there 
is not a large population of gas rich low luminosity or low surface brightness 
galaxies that has gone unnoticed by optical surveys. The infiuence of large scale 
structure on the determination of the HiMF from the Arecibo H I Strip Survey 
is tested by numerical experiments and was not found to affect the resulting 
HiMF significantly. The cosmological mass density of H I at the present time 
determined from the survey, ^lm{z = 0) = (2.0 ± 0.5) x 10^^h~^, is in good 
agreement with earlier estimates. We determine lower limits to the average 
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column densities (A'hi) of the galaxies detected in the survey and find that none 
of the galaxies have (A^hi) < lO^^'^cm"^, although there are no observational 
selection criteria against finding lower density systems. Eight percent of the 
signals detected in the original survey originated in groups of galaxies, whose 
signals chanced to coincide in frequency. 

Subject headings: galaxies: luminosity function, mass function - ISM - surveys 
- radio fines: galaxies — galaxies: ISM 
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1. Introduction 

The distribution function of neutral hydrogen masses among galaxies and intergalactic 
clouds (the H I mass function or HiMF), and more generally, the neutral hydrogen density 
in the nearby universe, Qm, are important inputs into models of cosmology and galaxy 
evolution. Different attempts have been made to construct an HiMF by using optically 
selected galaxies (Rao & Briggs 1993, hereafter RB, Solanes, Giovanelli & Haynes 1996). 
These studies are based on the assumption that H I is always associated with optically 
bright galaxies. A major concern is whether the HiMF is complete when it is computed for 
these galaxies. 

For example, the population of low surface brightness (LSB) galaxies might 
hypothetically constitute a substantial portion of the population of nearby extragalactic 
objects (Disney 1976, McGaugh 1996, Dalcanton, Spergel & Summers 1997, Sprayberry 
et al. 1997). The LSB population easily escapes detection optically and would not be 
included in the samples that are commonly used to evaluate the HiMF. This would be 
particularly problematic since LSB galaxies are generally found to be rich in neutral gas 
(Schombert et al. 1992, de Blok, McGaugh & van der Hulst 1996) and could therefore 
contribute substantially to the neutral gas content. Gas rich dwarf galaxies may also play 
an important part. For example Tyson & Scalo (1988) have argued that the majority of 
these dwarf galaxies remain undetected because only the small portion that is currently 
undergoing a rapid phase of star formation is presently observed with optical telescopes. A 
final possible population of gas rich systems that would escape inclusion in optical catalogs 
is a class of intergalactic H I clouds without stars. So far, only a few such systems have 
been discovered, and they are always found to be gravitationally bound to a galaxy or a 
group of galaxies (for example the Leo ring, Schneider 1989). 

Clearly, the HiMF and the H I content of the local universe should be measured 
directly, in such a way that they suffer no bias against gas rich galaxies or intergalactic H I 
clouds which are difficult to detect optically. This is possible by means of surveys in the 
H I line. Several of these surveys have been carried out over the last two decades. The 
majority were single dish observations using on/off techniques, since these surveys were 
done in conjunction with observations targeted on cataloged galaxies (Fisher & TuUy 1981b, 
Giovanelli & Haynes 1985, 1989). Some surveys were concentrated on groups of galaxies 
(Haynes & Roberts 1979, Lo & Sargent 1979, Fisher & TuUy 1981a, Hoffman, Lu & Salpeter 
1992) to specifically search for H I clouds in the vicinity of known galaxies. Other surveys 
were designed to find H I signals in voids (Krumm & Brosch 1984, Szomoru et al 1996), 
or to compare voids and superclusters (Weinberg et al. 1991, Szomoru et al. 1994). Since 
these surveys are not pointed at randomly chosen regions of sky, they may not provide fair 
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tests of the shape of the HiMF or an unbiased measure of the average H I density. A few 
truly bhnd surveys have been conducted, the first one in driftscan mode (Shostak 1977) 
and one by observing a series of pointings on hues of constant dechnation (Kerr & Henning 
1987, Henning 1992). It is worrisome that this latter survey could not reproduce the HiMF 
defined for optically selected galaxies, possibly because the survey was targeted toward 
large volumes of known voids. Another possibility is that the achieved survey sensitivity is 
not well understood, leading to an underestimation of the true number density of H I rich 
galaxies and intergalactic clouds. 

More recently, several surveys have been made using the Arecibo telescope (Sorar 
1994, Spitzak 1996, Schneider 1997). Surveys are in progress at Dwingeloo (DOGS) and at 
Parkes, where the survey will cover the entire Southern Sky (Staveley-Smith et al. 1996). 

In this paper we analyze the Arecibo H I Strip Survey (Sorar 1994), a blind survey 
for extragalactic H I covering ~ 65 square degrees of sky, out to a redshift of 7400 kms~^. 
The analysis of the survey results will concentrate specifically on the understanding of 
the achieved survey sensitivity and the vulnerability to large scale structure. We describe 
the details of the survey and the optical and 21cm follow-up observations in section 0. 
Section ^ gives a detailed analysis of the achieved survey sensitivity. We present the HiMF 
in section ^ the possible influence of large scale structure on the determination of the HiMF 
is examined by performing numerical experiments, and the cosmological mass density of 
H I at the present time is calculated in this section. In section ^ we compare our findings 
with previous estimates of the HiMF based on 21cm surveys and optically selected galaxy 
samples and discuss the implications of our results. Section ^ summarizes the results. The 
distances used in this paper are based on a Hubble constant Hq = lOOkms"^ Mpc~^. 

2. Observations 

2.1. Arecibo H I Strip Survey (AHiSS) 

The strip survey was carried out on the AreciboQ 305m Telescope in the period of 
August 1993 until February 1994. The survey was designed to take advantage of periods 
of construction when the telescope pointing was immobilized but the receiving systems 
were still operational. Therefore, the data were taken in driftscan mode and the telescope 
beam traced strips of constant declination. The same strips were retraced for as many as 



^The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated 
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30 days in order to obtain very sensitive observations that are capable of detecting H I of 
low surface density. The main survey was divided in two sessions: one survey covering 10.5 
hours of RA at 6 = 14°14', and the second survey covering 9.7 hours of RA at 5 = 23°09'. 
All observations were made at night. The limiting column density was lO^^cm"^ (at a 5a 
level) for gas filling the telescope beam, which subtends 3 kpc at 3 Mpc and 70 

kpc at 74 Mpc. The total sky coverage was approximately 65 square degrees, 
with a depth of cz = 7400 kms~^. The survey was capable of detecting H I masses of 
6 X lO^/i-^M© at 7 h-^ Mpc and 1.5 x lO^/i'^M© at the full depth of the survey, in the main 
beam which has a FWHM of 3 arcmin. The first sidelobe of the telescope beam pattern 
could detect 1.5 x lO^/i~^M0 at the full depth of the survey, which makes the sidelobes 
effective in detecting high H I mass galaxies over a 15 arcmin wide strip. The details of the 
AHiSS and the data reduction are described by Sorar (1994). A summary of the reduction 
path is described by Briggs et al. (1997). 

The survey yielded a total of 61 detections, of which approximately half could be 
associated with cataloged galaxies listed in the NASA Extragalactic Database (NED^. Five 
detections have no obvious counterparts on the Digitized Sky Survey (DSSQ) although they 
are more than 10° away from the Galactic plane. 

Fig. shows slice diagrams of the location of all H I selected objects compared with 
galaxies in the CfA catalog (Geller & Huchra 1989). Galaxies within a 10° strip centered on 
the declination of the strip are plotted. Different symbols are used to distinguish between 
cataloged (circles) and uncataloged (boxes) galaxies. The size of the symbols reflects the 
H I mass of the galaxies, in the sense that bigger symbols indicate larger H I masses. The 
Galactic plane intersects the strips at RA ^ 6^^ and RA ^ 19^^ as indicated by the dashed 
lines. The CfA catalog is obviously incomplete in these regions of the sky due to Galactic 
extinction, whereas the H I survey suffers no bias against detection in these regions. It is 
clear from this figure that H I selected galaxies generally follow the same structures as the 
optical galaxies. In particular, if we only consider those regions of the sky where the CfA 
catalog is complete, we find that more than 80% of the H I selected galaxies lie in regions 
where the average galaxies density is higher than the mean density. This is consistent with 



•^Thc NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space Administration. 

^Based on photographic data of the National Geographic Society - Palomar Observatory Sky Survey 
(NGS-POSS) obtained using the Oschin Telescope on Palomar Mountain. The NGS-POSS was funded by 
a grant from the National Geographic Society to the California Institute of Technology. The plates were 
processed into the present compressed digital form with their permission. The Digitized Sky Survey was 
produced at the Space Telescope Science Institute under US Government grant NAG W-2166. 
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the finding that LSB galaxies and gas-rich dwarfs he on structures dehneated by normal, 
high surface brightness galaxies (Bothun et al. 1986, Thuan, Gott & Schneider 1987, Eder 
et al. 1989, Thuan et al. 1991, Mo, McGaugh & Bothun 1994). Furthermore, none of the 
H I selected galaxies in the sections where the CfA is complete are found in regions where 
the galaxy density is less than one fifth of the cosmic mean. This is in agreement with the 
results of Szomoru et al. (1994) and Weinberg et al. (1991) that H I selected galaxies are 
not found in selected void fields. A more detailed analysis of the large scale distribution of 
the H I selected galaxies will be presented elsewhere. 





Fig. 1. — Slice diagrams of the location of all the H I selected objects from the Arecibo H I Strip 
Survey. Boxes indicate uncataloged galaxies, circles indicate cataloged galaxies. The size of the 

symbols reflects the H I mass of the galaxies. Also shown are galaxies from the CfA catalog (Geller 
& Huchra 1989) within a 10° strip centered on the declination of the surveys. The Galactic plane 
intersects the strips at RA« 6^ and RAss 19^ as indicated by the dashed lines. 
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2.2. 21cm Follow-up observations 

Follow-up 21cm synthesis observations on all signals found in this survey have been 
performed with the NRAOQ Very Large Array (VLA) in D-configuration. These follow-up 
observations are essential for three reasons: 

1) As stated before, the survey was capable of detecting signals as far as 7 arcmin from 
the center of the survey strip. Consequently, the coordinates yielded by the survey have 
uncertainties on the order of several arcminutes. The association with a cataloged galaxy or 
a galaxy on the DSS is therefore not always unambiguous, especially in those cases where 
there are several prominent galaxies visible within a 15 arcmin wide strip. 21cm observations 
with spatial resolution of ~ 1 arcmin, are sufficient to obtain unique identifications. 

2) Flux measurements from the AHiSS can be poor if a signal is detected at large distance 
from the center of the beam. In principle, a correction to the fiux could be made since 
the response function of the telescope is known with reasonable accuracy, but this is only 
possible if the positional accuracy is sufficient. Furthermore, measurements from the survey 
spectra might underestimate the integral flux if the source is more extended than the 
primary beam. Thus, a rough measure of the H I distribution increases confldence in the 
analysis. 

3) Some signals can be caused by pairs or small groups of galaxies, whose fine emission 
might stack up in the same channels. It is not obvious from the Arecibo survey spectra 
which signals are caused by more than one galaxy. In fact, this situation was found by the 
VLA observations to occur in flve cases. 

Short VLA observations (~ 20 min) of all 61 detected galaxies were performed during 
the D-conflguration sessions in May 1995 and September 1996. The signal of three systems 
fell below the detection limit of the snap-shot observation. These systems were re-observed 
in the D-conflguration during the second session, but with longer integration times (~ 3 
hours), resulting flnally in conflrmation of all 61 detections at levels consistent with the 
AHiSS sensitivity. The three weaker signals originated in galaxies whose declinations are at 
the center of the AHlSS. 

The VLA observations were performed using 63 channels over a 3.125 MHz bandwidth, 
corresponding to a velocity range of approximately 660 kms~^. On-line Hanning smoothing 
was used, resulting in a velocity resolution of ~ 10.5 kms^^. Phase calibrators were observed 
once for each source. Only seven different phase calibrators were used. Each time that the 
phase calibrator was changed, a primary flux calibrator (3C48 or 3C286) was observed in 



'"'The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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order to tie the flux scales together. Therefore, the primary cahbrators were not observed 
with the same correlator settings as each galaxy, but only with six correlator settings. The 
overhead time due to slewing of the telescopes and observations of primary calibrators is 
signiflcantly decreased this way. This technique provided adequate passband calibration for 
these short observations, where high dynamical range or accurate channel-to-channel flux 
density calibrations are not needed. There were no strong continuum sources in these flelds. 

Since the observing conditions were generally quite good, little editing was necessary 
to remove interference and bad baselines. Continuum was removed from the data in the 
Mf-plane, by making linear flts to the real and imaginary components for each visibility in 
the line-free channels and subtracting the appropriate values from all channels. The uv 
data were calibrated and transformed to datacubes using natural weighting. Using natural 
weighting rather than uniform weighting results in a slightly lower spatial resolution, but 
higher sensitivity. The resolution in the transformed datacubes is ~ 60" x 60". The flnal 
r.m.s noise was approximately 1.2 mJy/beam for each channel (corresponding to a minimal 
detectable column density of ~ 1.7 x lO^^cm"^ for the short integrations and 0.8 
mJy/beam (limiting column density ~ 9 x lO^'^cm"^ for the longer integrations. Since 
the observed galaxies were generally barely resolved, CLEANing of the data was not really 
necessary. Nonetheless, because the synthesized beam of these short observations has strong 
deviations from a Gaussian shape, we chose to clean the datacubes and restore them with 
a Gaussian beam in order to make an effective search for companion objects throughout the 
primary beam. 

Total H I masses were calculated using Mhi/M© = 236 d^/ SdV , where d is the 
distance to the source in Mpc, and S is the flux density in mJy over proflle width in 
kms"^ Total H I maps were constructed by adding the regions in the channel maps that 
contain line emission. Contour maps of H I emission and global proflles will be presented in 
a future paper. 

Fig. 1^ shows the H I masses as measured with the VLA plotted against the H I 
masses that are derived from the AHiSS spectra. Different symbols are used to distinguish 
between single detections (fllled circles) and multiple detections (open circles). Arecibo 
measurements of these multiple detections are of course always unreliable estimates of 
the real H I masses. This flgure clearly illustrates the fact that the H I masses derived 
from the Arecibo spectra generally underestimate the true H I masses, and that follow-up 
observations were essential to determine accurate measurements. In the analysis of the 
AHiSS results we will use H I masses calculated from the VLA observations, except for 
galaxies with integrated VLA fluxes less than 1.0 Jykms"^. Galaxies with these low 
fluxes are only found at small distances from the center of the Arecibo beam, where the 
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normalized response function is close to unity. Arecibo measurements for these galaxies are 
therefore reliable estimates of their real fluxes. 

Table |l| summarizes the global parameters of the detected galaxies which are 
derived from the VLA observations. The following information is contained: Column 
1: Identification number of the Arecibo detection. Indices indicate multiple detections. 
Column 2: Name of galaxy if already cataloged. Columns 3 and 4: B1950 coordinates. 
Column 5: Logarithm of H I mass. Column 6: Heliocentric velocity, calculated by taking 
the mean of the velocities at 20% of peak flux density. Column 7: Declination offset from 
center of the survey strip. Column 8: Identification code for multiple detections. 
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Fig. 2. — Comparison of total H I masses as measured with the VLA and H I masses derived from 
the Arecibo H I Strip Survey spectra. The dashed line is the line of equality. Different symbols 
distinguish between singe detections (filled circles) and multiple detections (open circles). 



2.3. Optical follow-up observations 



We have also started a program of optical i?-band imaging of the H I selected galaxies 
with the 2.5m Isaac Newton Telescope^ (INT) on La Palma. During four observing runs 
in the period between October 1995 and March 1997 we have been able to obtain images 
of all target galaxies. The first noteworthy result of the follow-up optical observations is 
that all of the H I sources more than 10° away from the Galactic plane are found to be 
associated with galaxies in the 5-band images. We have seen no indication that any of the 
sources detected in the Arecibo survey are anything other than 'ordinary' galaxies having 
both gas and stars. Put another way, we have so far failed to find any H I sources that are 
pure H I clouds without stars. The analysis of the optical observations will be presented in 
a subsequent paper. 



3. The survey sensitivity 

The improved coordinates and H I fiuxes derived from the VLA observations help to 
verify the completeness limit of the survey. The detectability is determined in the following 
way: 

The H I mass in a detected signal can be expressed as Mhi/M© = 236 x S AV(P, 
where d is the distance to the source in Mpc, and S is the fiux density in mJy that here 
is considered to be constant over a rectangular profile of width AV^ km/s. The sensitivity 
of an observation is optimal when the spectra are smoothed to the velocity width of 
the source. With optimal smoothing the noise becomes a{AV) = ao^J AVo/ AV , where 
AVo is the spectral resolution of the receiving system and ctq is the noise level in the 
unsmoothed spectra. The limiting fiux density Sc for a 5a detection is then given by 
S'c(AV") = 5a{AV) = 5ao\fAVojAV. For the AHiSS, the average noise level after coadding 
spectra taken at different days was 0.75 mJy for a velocity resolution of 16 km/s. 

The normalized response function of the survey telescope, I{0), describes the relative 
response to a source which is detected at an offset 6 from the center of the survey strip. In 
other words, I{6) is the integral of the flux density sensed by the telescope as a source makes 
a cut through the beam pattern, missing the center of the beam by angle 6, normalized in 
such a way that /(O) = 1. For the Arecibo telescope this function falls off to ~ 0.125 at 
6 = 3.25' and reaches a second maximum of ~ 0.2 at ^ = 4.8' due to the high sidelobe level. 



^The Isaac Newton Telescope is operated by the Royal Greenwich Observatory in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto Astrofisica de Canarias. 
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We define a correction factor Cj-{9) that accounts for the shape of the response function by 
c^{6) = 1/I{6). The hmiting fiux density of a galaxy at an offset 6 from the center of the 
survey strip can then be expressed as ScC^{9). In general, an H I source should be detected 
by the survey if its flux density exceeds this limiting flux density, that is 5* > 5'cCr(^). This 
can be rewritten as 



where we deflne D as the detectability. The factor Cf represents the normalized feed gain 
of the telescope, which is a function of frequency. Cf can be approximated by an analytical 
expression of the form Cf ~ 1 — ([/ — /o]/w)^, where / is the frequency in MHz, /o is the 
center of the survey band and w is a parameter which determines the width of the band. 
During the observations at 5 = 14° 14' the shape of the normalized feed gain remained 
unchanged and could be satisfactorily flt by /o = 1395 MHz and w = 40.5 MHz. During the 
S = 23°09' observations, the gain was retuned a few times to the settings: /o = 1395 MHz 
and w = 37.5 MHz and /o = 1410 MHz and w = 33.8 MHz. The Cf dependence of D leads 
to a 'distance dependence' for the flux density sensitivity. 

The sensitivity can now be verifled by plotting the detectability D of the objects 
against declination offset from the center of the survey strip. This is shown in Fig. ^ All 
signals that were detected are shown in this plot. The horizontal errorbars are the result 
of a combination of positional accuracy in the VLA maps and the spatial extent of the 
galaxies in the direction orthogonal to the survey strip. The signals corresponding to pair 
or group detections are indicated by open circles and letters A to E are used to identify 
these related or confused signals in Table |I|. Filled circles mark all single object detections. 
The solid line represents Cr(6'), the limit to the detectability. In principle, all fllled circles 
should lie above the solid line, the area below this line is 'undetectable'. The line is a 
satisfactory limit to the data points, especially if we consider the naive character of Eq.^. 
That is, this equation assumes that the detected proflles are symmetric and featureless. 
Since the proflles are generally heavily smoothed, this is a reasonable assumption, but 
lop-sided proflles or strong double horned proflles might exceed the detection limit while 
they are formally 'undetectable' according to Eq.|l[ 

The limiting depth d^, the maximum distance to which the object could be placed and 
still remain within the sample, can now be expressed as a function of Mm, 9 and AV: 




(1) 



4(MHi,e,Ay) 



Mhi/(^) 



(2) 



\ 236S'c(A\/) AV 



The variation of feed gain with frequency imposes a minor correction to the limiting 
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Fig. 3. — Detectability, D, of the objects versus the offset from the center of the survey strip. 
Multiple detections are indicated by open symbols and coded to indicate which galaxies originate 
from the same Arecibo detection. The solid curve shows the detection limit. Objects below this 
line should be undetectable if their H I profile widths would be rectangular. 



depth dc which can be compensated by solving for dc in 

S{dc)C{{dc) = Sc, (3) 

that is, the true flux density S{dc) of a galaxy located at the distance limit dc multiplied by 
the feed gain at a frequency corresponding to that distance has to be equal to the limiting 
flux density Sc obtained at the center frequency, where the feed gain is optimized and the 
nominal noise level of the system is normally calculated. 

The next step is to calculate the effective search volume of the survey. The volume dV 
of a slice of the total survey volume at declination 6 and with length / (the total length of 
the strips in radians, / = ARA27rcos5/24'^) and width d9 is given by 

MM,,e.Av) = I if^'^:'''^^''"- i!t " (4) 

where dsw is the limiting depth of the survey imposed by the bandwidth of the receiving 
system. The total survey volume that is sensitive to a galaxy with H I mass Mhi and 
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velocity spread AV^ can then be calculated by taking the integral over 9. For galaxies in the 
flux limited regime, for which < c?bw for each 6, this integral simply becomes: 

For galaxies in the bandwidth limited regime, dc > dsw for ^ smaller than a certain critical 
value 6'bw- The integral must now be split up in separate parts for the flux limited regimes 
and the band width limited regime: 

1 /•-6'bw „ 1 rdBw „ 1 r°° n 

v(Mhi,av) = - / dlide + - divide + - dlide 

where we made use of the symmetry of the beam shape in the last step. 

This effective search volume is still dependent on two variables: Mhi and AV. The 
most convenient parameterization of V for computing an HiMF is to express V as a function 
of Mhi only. This can be achieved by adopting a relation between Mm and AV^. Such a 
relation is known to exist since optical luminosity L is related to AV via the TuUy-Fisher 
relation (TuUy &; Fisher 1977), and L is related to Mhi as Mhi oc (see Briggs 1990). 
Briggs & Rao (1993, hereafter BR) determined the M-m-AV relation empirically by plotting 
AV against logMni for 1139 optically selected galaxies from the catalog by Fisher & Tully 
(1981b). A fit to these points gives AV = 0.16 Mh{^ Recently, Salpeter & Hoffman (1996) 
analyzed H I observations of 70 dwarf galaxies and find a similar trend: AV oc M^^- 
This relation is therefore valid over a wide range in H I mass. Our data are also in good 
agreement with this relation. Note that the velocity widths in these relations are not the 
inclination corrected maximum rotational velocities, but just the observed velocity spreads. 

The effective survey volume V(Mhi) can now be calculated by substituting the relation 
between AV and Mhi found by BR and Sc{AV) into Eq.^ and Eq.^. This volume as a 
function of Mhi is shown by the solid line in Fig. ^. In the flux limited regime, for H I 
masses < 10^'^, the search volume is V oc Mhi^. (A proportionality often used in the 
literature is V oc M^{^ [e.g., Henning 1995, Schneider 1997]. This power 3/2 arises if the 
dependence of AV on Mhi is discarded, effectively assuming that all galaxies have the same 
profile width.) In the high mass region (Mhi > 10^'^), the limiting depths are no longer 
determined by the detectability of the signals, but simply by the bandwidth of the receiving 
system. Therefore, the effective survey volume for high mass systems is not dependent on 
H I mass or AV. The total survey volume in this regime, is ~ 3000/i-3Mpc^ For lower 
H I masses the volume decreases rapidly and is only ~ 1.0/;.~^Mpc'^ for Mhi = lO^M©. The 
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dashed and the dotted line show the effective search volume corresponding to the main 
beam and the sidelobes, respectively. This figure clearly shows that the sidelobes do not 
add much volume in the low mass range, but are very effective in finding large H I masses. 




Fig. 4. — The effective search volume vs. H I mass. The solid line is the total search volume 
of the Arecibo H I Strip Survey, the dashed and the dotted line show the effective search volume 
corresponding to the main beam and the sidelobes, respectively. The arrows indicate the survey 
volumes of two H I surveys which are comparable in size to the AHiSS: Henning (1995) and Spitzak 
(1996). 

Since the observed velocity width is related to inclination i and maximum rotational 
velocity Vmax as AV — 2VinaxSini, the hmiting depth is dc oc (AV)~^^'^ oc (sini)~^/^. The 
limiting depth is therefore weakly dependent on inclination, in the sense that low inclined 
galaxies can be detected to slightly larger distances. Consequently, an H I survey would in 
principle preferentially select face-on galaxies. In practice, this effect is negligible, since the 
expectation for ((sini)~^/^) of a randomly oriented sample is close to unity: 

Jo sm I di 



Only the flux limited regime would be hampered by the inclination effect, but most galaxies 
in this regime are low mass, dwarf galaxies. The velocity spreads are small for these 
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galaxies, and turbulent motions play an important part at establishing the profile width 
(cf., Lo, Sargent & Young 1993). 

A potential hazard in radio surveys is the influence of radio frequency interference 
(RFI). RFI can in principle cause false negatives (miss a signiflcant signal) in the sample if 
it affects the spectrum exactly at the frequency of a source, or when a narrow-lined source 
is mistaken for RFI. However, the driftscan method that has been used for the AHiSS has 
proven to provide a very good stability for RFI signals (see Briggs et al. 1997). Repetitive 
coverage of the same regions of sky makes the survey immune to RFI and unstable baselines. 
This has been demonstrated by the fact that all signals identifled in the Arecibo survey 
have been conflrmed by the VLA follow-up. 



4. The H I mass function 

The H I mass function (HiMF) is deflned analogously to optical luminosity functions. 
The HiMF Q{Mm/M^i)d{Miii/M^i) gives the total number of galaxies or intergalactic 
clouds per Mpc^ in the mass interval (i(MHi/Myj) centered on Mhi/M^j. Here, we flnd 
it convenient in our analysis and flgures to plot the HiMF as the number of galaxies or 
intergalactic clouds per decade in mass. In order to parameterize the shape of the HiMF, 
we adopt the conventional Schechter (1976) function, 

e(^M(^) = r(^)-"exp-(^M^*!i), (8) 

with free parameters a, the slope of the low- mass end, Mj^j the characteristic mass that 
the deflnes the kink in the function and 6*, a normalization factor. This function must be 
integrated over decade bins for comparison with the binned data in the flgures. 



4.1. Methods 

The classical method (see e.g., Christensen 1975, Schechter 1976) for determining 
luminosity functions is based on the assumption that galaxies are distributed in a uniform 
manner. The luminosity function is determined by dividing the number of galaxies in a bin 
centered on M by V(M), the effective search volume for that particular M. This method is 
easily applicable to the AHlSS. The search volume can be evaluated with Eq.^ and Eq.^, 
using a statistical relation between velocity width and H I mass. The advantages of this 
method are that it is nonparametric and that it is automatically normalized. The important 
disadvantage is that it assumes homogeneity and its use might lead to errors in 9(Mhi) if 
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density fluctuations due to large scale structure occur on distance scales comparable to, or 
greater than the depth at which Mm can be detected. 

A slightly modified form of the 'classical' method is the I](l/Vmax) method, first used 
by Schmidt (1968), which is also applicable to H I surveys. Instead of calculating a mean 
survey volume for each mass bin, this method consists of summing the reciprocals of the 
volumes corresponding to the maximum distances to which the objects could be placed 
and still remain within the sample. The values of Vmax can be calculated directly using 
Eq.|^ and Eq.||, using now the measured velocity spread instead of the statistical value. 
The two procedures give similar results when the number of galaxies per bin is large. For 
the less densely populated bins the two procedures can give different results because Vmax 
of a particular galaxy can strongly deviate from the average V of the mass bin it falls in. 
In section [4.2| the HiMF is calculated using Schmidt's method. Like the classical method, 
Schmidt's method is also vulnerable to errors caused by large scale structure. The possible 
effects of large scale structure are discussed in section ^ISj 



4.2. Results 

Fig. 1^ shows the principal results of this analysis. The lower panel shows the observed 
distribution of H I masses binned per half decade, with errorbars given by Poisson statistics. 
This histogram shows that the survey has detected galaxies with H I masses in the range 
from 10^ to IO^^Mq and therefore enables us to determine the HiMF over three orders of 
magnitudes in H I mass. 

The solid dots in the upper panel of Fig. |^ show the HiMF determined by the 
X](l/Vmax) method. The errorbars are given by Poisson statistics. Also drawn in this figure 
are analytical curves given by the Schechter function of Eq A satisfactory fit to the points 
is obtained with a = 1.20, 9* = 0.014 Mpc'^ and log(Mj^j/Mo) = 9.55. Mass functions 
with faint end slopes of 1.10 and 1.30 are shown to indicate the uncertainty in the value of 
a. We note that in the present analysis the parameterization of the HiMF in the form of 
a Schechter function is only used to enable comparison with other H I survey results and 
results based on the optically selected galaxy population. The I](l/Vmax) method recovers 
the shape and amplitude of the HiMF simultaneously without using a Schechter function 
(or any other parameterization) as an assumption about the intrinsic shape of the HiMF. 

The observational limits to the determination of the HiMF are also illustrated in Fig. ^ 
by the thin line. This line represents the sensitivity function of the survey to objects of H I 
mass Mhi, defined as 0(Mhi) = 1/V(Mhi), where V(Mhi) is the effective search volume. 
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Fig. 5. — Lower panel: The distribution of H I masses of the detected galaxies from VLA follow-up 
measurements. The errorbars are given by Poisson statistics. Upper panel: The thin line is the 
sensitivity of the survey defined by </> = 1/V, where V is the effective search volume. In the region 
IO^Mq < Mhi < IO^^Mq this function defines an upper limit to the space density of intergalactic 
H I clouds without stars. The measured H I mass function per decade is shown by the points. The 
fat line is a Schechter luminosity function with parameters as given in the upper right corner. Also 
Schechter functions with a = 1.1 and a = 1.3 are shown. The arrows show upper limits to the 
space density of H I rich galaxies or intergalactic H I clouds. The two arrows on the right are from 
a complementary survey with the Arecibo telescope over the range 19,000 to 28,000 kms~^. 

In the ranere IO^Mq < Mhi < 10^°Mq this function defines an upper limit to the space 
density of intergalactic H I clouds without stars. It also shows that this survey is not 
capable of measuring the HiMF directly in the regions Mri < lO'^'M© and Mm > lO^^M© 
if the extrapolation of the analytical form for ©(Mhi) holds. The sensitivity function does 
allow us to define upper limits in these ranges as indicated by arrows. The upper limit at 
Mhi = 10^°'^'''Mq is from the complementary survey with the Arecibo telescope over the 
range 19,000 to 28,000 kms~^, only sensitive to high mass galaxies (See Sorar (1994) for 
details) . 



- 18 - 



4.3. Influence of large scale structure 

The slice diagrams in Fig |1] show that the survey strips sample a wide range of large 
scale structure, as the combined RA range extends nearly 2/3 of the way around the sky. 
On the other hand, the limited depth of the survey for small H I masses might cause 
the low mass end of the HiMF to be affected by local density fluctuations. Therefore, 
a major concern in the determination of the mass function following Schmidt's method, 
is whether homogeneity is a fair assumption. For instance, a local density enhancement 
would overestimate the number of low mass galaxies or clouds and would give rise to 
an overestimate of the faint end slope of the mass function. The biases in the shape 
of the mass function due to large scale structure can be avoided by making use of a 
maximum likelihood estimator (cf., Efstathiou, Ellis & Peterson 1988, Saunders et al. 
1990). If the spatial density distribution of galaxies and intergalactic clouds is given by 
N{Miii,r), then the overall density p{r) and the mass function 6(Mhi) can be separated as 
N{Mi{i,f) = p(r*)B(MHi). Maximum likelihood estimators can then be used to solve either 
for O(Mhi) without knowledge of p(r) or solve for p(r) without knowledge of 0(Mhi) . 

Although the maximum likelihood method is a powerful procedure for determining 
optical luminosity functions (Marzke, Geller & Huchra 1994, Lin et al. 1996, Saunders et al. 
1990), we chose not to apply it to the AHlSS. The most important problem is that we have 
to deal with small number statistics. Especially the faint end slope of the mass function is 
deflned by very few galaxies or clouds per bin. Some experimentation with application of 
the algorithm to small samples such as ours showed that maximum likelihood methods can 
produce erratic results is these situations. Another complication arises because maximum 
likelihood methods assume that the shape of the H I mass function is independent of the 
space density p(r). In practice, high mass galaxies will have a higher statistical weight 
in the determination of p{r) since they are simply more numerous in our sample. It is 
questionable whether this space density deflned by the high mass galaxies is a fair estimate 
for the galaxies at the faint end side of the mass function, in the region that is dominated by 
dwarf and LSB galaxies. Although these galaxies are found to follow the same general large 
scale structures as the normal HSB galaxies, they preferentially avoid the highest density 
regions (Mo et al. 1994, Taylor 1996). In principle, the spatial density p could be separated 
into different functions for individual morphological types or different ranges in H I mass, 
but the limited number of galaxies in our sample does not allow this differentiation. 

In order to investigate the effect of large scale structure we have performed numerical 
experiments. These tests consisted of randomly placing artiflcial galaxies in a volume with 
weights given by the HiMF given in Eq.^ A range of HiMF parameters were investigated. 
The results are illustrated using the derived HiMF parameters from Fig^ 6 = 0.014 Mpc~^, 
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log(Mj^j/M0) = 9.55 and a = 1.20. The galaxies have random inchnation, and rotational 
velocity related to H I mass as determined by BR. Galaxies were selected from these 
volumes in the same manner as the Arecibo H I Strip Survey selects galaxies from the sky. 




Fig. 6. — Numerical tests of the influence of large scale structure on the determination of the HiMF. 
The left panels show the imposed normalized density fluctuations in the artificial catalogs. The 
histograms show the redshift distributions of the simulated data as solid lines, with la uncertainties 
as thin lines. The measured redshift distribution of our survey sample is represented by solid circles 
and a dashed line. The right panels show the J2{^/^ma.x) HiMF of the simulated data as open 
squares, with la err or bars. The input HiMF is indicated by the solid curves. 

Fig. D shows the results of three experiments, labeled as A, B and C. The presented 
results are each the average of 100 independent simulations. The top panels show the 
assumed normalized density fluctuations p. The central panels show the redshift distribution 
of the simulated data as thick-lined histograms, thin solid lines show la uncertainties. 
Overlaid on these histograms is the measured redshift distribution of our survey sample, 
indicated by solid circles and dashed histograms. The bottom panels show the averaged 
HiMFs of the simulated catalogs as open squares, with la errorbars computed as the 
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standard deviations of the 100 simulations. The Yl 1 /Vmax niethod has been used. 

In case A a homogeneous density distribution is assumed. It is not surprising to see 
that the input HiMF and the measured HiMF agree with high accuracy. The 
method should give reliable results in this situation. It is however apparent that the survey 
observations have several distance bins that are overdense or underdense compared to the 

uncertainty band of the 'uniform density' simulation. The comparison implies that we may 
be observing underdensities at D ^ 30 Mpc and D ^ 70 Mpc and an overdensity at D ^ 50 
Mpc. To test this hypothesis we constructed catalogs with a density fluctuation as indicated 
by the top panel in case B. The precise functional form of this fluctuation is not important, 
the condition is that is should produce consecutive under and overdensities at D ~ 30, 50 
and 70 Mpc. The central panel shows that the imposed density fluctuation reproduces the 
observed redshift distribution satisfactorily. It is surprising however that if the Z^l/Vmax 
method is applied to these simulated data the resulting HiMF is still indistinguishable from 
the input HiMF. In other words, the J2 1/Vmax method appears to be a robust method, and 
not very sensitive to the effects of large scale structure. This is even true for case C, where 
a more severe density fluctuation is imposed, and where the resulting HiMF is still in good 
agreement with the input HiMF. We can conclude from these simulations that although we 
see a hint of the effects of large scale structure in our data, the observed deviations from 
uniformity have no significant influence on our determination of the HiMF. 



5. Discussion 

5.1. Previous estimates of the HiMFs from H I surveys 

The first H I surveys could only be used to set upper hmits to the space density of 
intergalactic H I clouds without stars, and did not yield enough detections to allow the 
determination of the shape of the HiMF. Shostak (1977) was to first to define these limits 
by means of blind driftscan emission and absorption searches in the 21cm line with the 
NRAO 91m telescope. Due to the relatively poor sensitivity (12-40 mJy) and the small 
total effective search volume (~ 300/i^'^Mpc^), these limits were not very strict. In the 
region Mhi < lO^'^M© the limits set by the Arecibo H I Strip Survey are at least two orders 
of magnitude lower (99% confidence) than those set by Shostak. Haynes & Roberts (1979), 
Lo & Sargent (1979) and Fisher & TuUy (1981a) search for invisible galaxies in groups of 
galaxies. Haynes & Roberts conclude from observations in the direction of the Sculptor 
group that intergalactic H I clouds with Mhi > IO^Mq do not exist. Lo & Sargent find four 
previously uncataloged LSB dwarf galaxies but their upper hmits to the space density of 
unseen objects do not improve Shostak's. Fisher & TuUy used the NRAO 91m telescope to 
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search for invisible galaxies in the M81 group. Their null result allowed them to push the 
upper limits to the space density 0.6 dex lower those set by Shostak. Also Krumm & Brosch 
(1984) find no H I sources in their driftscan searches of void regions and are only able to 
define upper limits. These limits are not very strict since their survey is only sensitive 
in the redshift range above 5300 kms~^, and is consequently only capable of finding H I 
masses > IO^^Mq. 

A series of papers by Kerr & Henning (1987) and Henning (1992, 1995) describes blind 
surveys which consisted of observing a series of pointings along lines of constant declination 
in the zone of avoidance and at high galactic latitude. These surveys yielded 39 detections 
(of which half were previously unknown) and were the first to put serious constraints 
on the shape of the HiMF. For comparison, the HiMF calculated by Henning (1995) is 
reproduced in Fig. |^ as solid squares, together with the AHiSS HiMF represented as a 
solid line. Henning's points are significantly lower than our measured curve. According to 
Henning's calculations objects with lO'^'^Mo < Mhi < lO^'^M© are deficient by a one order 
of magnitude compared to the AHiSS HiMF. The presence of the Local Void in Henning's 
survey could explain a part of this discrepancy, but even if the void is omitted the points 
are still beneath our HiMF. The most reasonable explanation for the discrepancy is that the 
sensitivity of Henning's survey is not well understood, leading to a serious underestimate of 
the HiMF. 

The most recent determination of the HiMF is the one by Schneider (1997) who used 
both the uncorrected Arecibo data from Sorar's (1994) survey and the data from Spitzak's 
(1996) survey with the Arecibo telescope. Schneider assessed the completeness limits of 
the surveys with the use of V/Vmax tests. He calculated the average value of V/Vmax and 
scaled that to the correct value of 0.5 by increasing the limiting fiux density Sc- Schneider 
noted the vulnerability of this method to the infiuence of large scale structure. Despite this 
different approach, the resulting HiMF is in good agreement with our estimate, as can be 
seen in Fig. |^. Only the low mass end deviates significantly. The statistics in the bins with 
Mui < 7.25 are poor, and Schneider could not rule out a slope as high as a = 1.7. The 
differences in the slopes can be partly explained by the fact that uncorrected H I masses 
are used in Schneider's estimate of the HiMF from Sorar's data. The effect that the flux is 
underestimated when the Arecibo beam misses a galaxy by a certain angle 6 was not taken 
into account. The VLA H I masses that are used in our HiMF are on average 50% higher, 
leading to a shift of galaxies to higher H I mass bins. As an alternate approach, Schneider 
also applied a maximum likelihood fit to all the signals brighter that 20cr. The resulting 
HiMF had a faint end slope a = 1.32. This result is also drawn in Fig. |^ and is in very 
good agreement with our function over the whole range in H I mass. 
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Fig. 7. — Comparison of H I mass functions of different surveys. The fat solid line shows the HiMF 
presented in this paper, using VLA observations of galaxies found in the Arecibo H I Strip Survey. 
Briggs k, Rao's (1993) HiMF based on the Fisher & Tully catalog is indicated by dashed curve. The 
solid squares are from Hcnning (1995). The open symbols result from Schneider's (1997) V/Vmax 
method applied to the data of Spitzak (1995) (circles) and the uncorrected Arecibo data from Sorar 
(1994) (squares). Schneider's (1997) maximum likelihood solution is indicated by the dotted curve. 
The dash-dot line shows a faint end slope of a = 1.7. 

In addition to making comparisons with surveys that are equally blind as the AHiSS, 
it is also interesting to compare our results with surveys that cover comparable volumes 
in underdense and overdense regions. It is noteworthy that the HiMF determined from 
Hi observations in clusters (McMahon 1993). as well as that determined in cosmic voids 
(Szomoru et al. 1996) are both well fit by a Schechter function, with approximately the 
same shape as as our HiMF. The scaling 9*, however, differs by a factor of ten or more 
between these different regions. 

Since the AHiSS crosses the Galactic plane twice, and consequently a large fraction 
of the total survey volume lies in regions which are not complete in optical catalogs, it is 
difficult to calculate the exact void filling factor of our survey. If we limit the calculation 
to the regions where |6| > 30°, we estimate that approximately 50% of the survey volume 
consists of regions where the average galaxy density is less than one third of the cosmic 
mean, and approximately 10% of the volume consists of structures where the galaxy density 
takes on more than three times the cosmic mean. The same percentages are obtained if 
random slices through the Universe are taken. If voids would be omitted from the survey, 
the scaling 9* of the HiMF would increase by a factor of two. 
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5.2. Comparison with HiMFs bases on optically selected galaxies 

The construction of an HiMFs from a sample of optically selected galaxies can either 
be done statistically, by studying optical luminosity functions and the dependence of H I 
mass on on optical luminosity, or directly by using 21cm data of optically selected galaxies. 
The first method was used by RB, who reviewed the literature that describes luminosity 
functions and H I richness for individual morphological types. Combining these data 
enabled them to construct an HiMF over three orders of magnitude in H I mass. Their 
result is in good agreement with our estimate of the HiMF in the range of H I mass where 
our survey is sensitive. 

The second method is discussed by BR who analyzed H I observations drawn from the 
catalogs of Fisher & TuUy (1981b) and Hoffman et al. (1989), and recently by Solanes et 
al. (1996) who use 21cm data of an optical magnitude-limited sample of galaxies in the 
direction of the Pisces-Perseus supercluster. Both authors arrive at the same values for 
the normalization and Mhj, but their values of the faint end slope differ significantly. BR 
find a ~ 1.25 while Solanes et al. find a declining slope. The Solanes et al. result may not 
be relevant to the discussion here since their sample excluded all dwarf irregular galaxies 
and contained no galaxies with Mhi < 10®'^ ~ ^Mj^j. The HiMF determined by BR is 
reasonably consistent with with our estimate of the HiMF as shown in Fig. |^. 

5.3. Implications: A new H I selected galaxy population? 

The comparison between our result and the HiMFs based on optically selected galaxy 
samples provides a direct test of the existence of a new population of gas rich galaxies that 
has gone unnoticed by optical surveys. In a way, this comparison appraises the completeness 
of optically selected catalogs for gas rich galaxies. Any contribution of uncataloged gas 
rich dwarf or gas rich LSB galaxies would yield an difference between HiMFs computed 
from optically selected and Hl-selected galaxy catalogs. Fig |^ shows that the HiMF derived 
from the results of the AHiSS is in very good agreement with previous estimates based on 
optically selected galaxy catalogs. This implies that the optically selected samples that 
have been used to evaluate the HiMF are remarkably complete. There is no evidence for a 
large number of neutral gas rich objects that have escaped inclusion in these catalogs. To 
the extent that optically selected catalogs are incomplete for LSB galaxies, the excluded 
galaxies must be predominantly gas poor, consistent with the finding by Sprayberry et al. 
(1997) that these excluded LSB galaxies are predominantly of low optical luminosity. 

Several authors have speculated about the existence of a large class of gas rich dwarf 
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galaxies (e.g., Dekel & Silk 1986, Tyson & Scalo 1988). These galaxies would yield a steep 
rise in the HiMF below Mhi ~ lO^M©. We find no observational support for the existence 
of this class of galaxies. Based on the counting statistics in the mass bins 10^'^ to lO^M© 
we can exclude a faint end slope of a = 1.7 or steeper with a 99% confidence level. 

5.4. Neutral gas density 

Our knowledge of Qg at high redshift is determined by the statistics of Damped Ly-a 
systems, seen in absorption against background quasars (e.g., Lanzetta, Wolfe & Turnshek 
1995, Storrie-Lombardi, Irwin & McMahon 1996). The picture emerging from these studies 
is that Qg reaches a maximum at 2; ~ 3 and has been declining since then. At low redshifts 
different effects complicate the determination of Qg. Firstly, at redshifts z < 1.6, the Ly-a 
line is not redshifted to optical wavelengths and has to be observed from space (Lanzetta 
et al. 1995, Rao, Turnshek & Briggs 1995). Secondly, the evaluation of fig may depend on 
the selection effects in the sample of quasars that has been used. Especially at these low 
redshifts, it is very difficult to compile a unbiased sample of quasars since QSOs would need 
to be observed within the optical images of galaxies. At the one hand, gravitational lensing 
can bring faint quasars into the sample which should otherwise be below the detection limit 
(Smette, Claeskens & Surdej 1997). But on the other hand, dust in Ly-a systems might 
obscure the background objects to a level where they are undetected (Fall & Pei 1993). 

The Arecibo H I Strip Survey can be used to evaluate the gas density differently. The 
mentioned problems are circumvented since fig is not measured by using the Ly-a line, but 
the 21cm line in emission. 

The space density of H I contained in objects of different H I masses is plotted in 
Fig. p. The solid line indicates an analytical expression of this function, determined from 
the product Mhi0(Mhi). The thin line represents the sensitivity limits and the arrows 
mark upper limits determined analogously to the upper limits in Fig ^. 

The integral H I mass density at the present epoch can be determined by taking 
the integral over the solid line in Fig. ^ This yields: pni = r(2 — a)6*M^i, where F 
is the Euler gamma function. Using the best fit Schechter parameters, we derive that 
Pjjj(^ = 0) = 5.8 X lO'^ h Mq Mpc"^ or 3.9 x lO'^^ h g cm'^, with a statistical error of 20%. 
A summation over all survey galaxies pni = J2 Afni/Vmax yields a shghtly smaller value: 
Pmiz = 0) = 5.4 X 10'^ HMq Mpc~^ since this calculation does not include the contribution 
to Phi of the density function below Mhi = IO^Mq or above Mhi = lO^^M©. The cosmological 
mass density of H I at 2; = is flmi^ = 0) = (2.1 ± 0.4) x 10"^/i~^ The total cosmological 
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Fig. 8. — Space density of H I mass contained in objects of different masses per decade. Tlie fat 
line shows the converted analytical H i mass function calculated by multiplying G x Mm, where Q 
is the HiMF plotted in Fig. ^. The dashed regions of the curve indicate that the sample contains 
no galaxies with Mhi < 10 ''Mq or Mhi > IO^'^Mq. The thin line indicates the sensitivity of the 
survey. The arrows mark upper limits to the space density of H i mass. The three arrows on the 
right are from a complementary survey with the Arecibo telescope over the range 19,000 to 28,000 
kms""*^. Note that galaxies with IO^Mq < Mhi < IO^'^'^Mq contribute most to the integral H i 
density. 

mass density of neutral gas at the present epoch is flg{z = 0) = (2.7 ± 0.5) x 10~'^h~^, 
assuming that the mass percentage of He I is 25% of the total gas density. 

This value agrees surprisingly well with earlier estimates by RB, who find 
Pmiz = 0) = (4.8 ± 1.1) X IO'^/iMq Mpc~^ by using optically selected galaxies. Fall 
& Pei (1993) arrive at approximately the same value of Qm by simply computing the 
average Mm/ Lb of the Huchtmeier & Richter (1989) catalog and multiply that by the 
mean luminosity density in the local universe as estimated by Efstathiou et al. (1988). The 
agreement between the cosmo logical mass density based on optically selected galaxies and 
that based on H I selected galaxies illustrates once more that there is not much neutral gas 
hidden in objects like LSB galaxies, dwarfs or intergalactic H I clouds that are missed by 
optical surveys. 

The contribution of dwarf galaxies to the H I content of the nearby universe is modest, 
galaxies with H I masses < IO^Mq make up only ~ 17% of the integral H I density. The 
density function in Fig. ^ clearly illustrates that the integral H I mass density is dominated 
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by high mass galaxies with H I masses in the range IO^Mq < Mhi < lO^^M©, which are L* 
galaxies. At Mhi ~ IO^'^Mq the density function drops off sharply, indicating that Malin 
1 type galaxies make no significant contribution to Ohi- This sharp cutoff was already 
noted by Bothun (f985). Much stronger upper limits to the contribution of galaxies with 
Mhi > 1O^°M0 wiU be set in the near future by the results of the Parkes Multibeam Survey 
(Staveley-Smith et al. f996). 

The estimate of the integral gas density from the AHiSS is a robust result. Two 
effects cause a relatively low uncertainty in the determination of Ohi- Firstly, the peak 
in the H I gas density function is conveniently caused by the galaxies that dominate the 
counting statistics. Galaxies in the lower mass bins, where the Poisson errors are large, 
contribute not much to the total density and therefore also not much to the uncertainty in 
Qhi- Secondly, the effective search volume in the mass region that dominates the density 
function is mostly band width limited. Uncertainties in Mhi and velocity width do not 
influence the search volume in this regime. 



5.5. What could be missed? 

The question arises which gas rich systems could be missed by the survey and could 
still contribute significantly to the integral gas mass density in the local universe. The only 
possible candidates that might escape detection are extremely low gas density systems with 
H I column densities below lO^^cm"^, the detection limit of the AHlSS. 

It has been shown that the gas density of spiral galaxies is correlated with the optical 
surface brightness in such a way that lower optical surface brightness implies lower gas 
densities. De Blok et al. (1996) observed a sample of LSB galaxies and showed that 
the neutral gas densities are generally a factor of two lower than those of late type high 
surface brightness galaxies. Furthermore, van der Hulst et al. (1993) has shown that the 
gas densities of LSB galaxies arc often just below or equal to the critical density for star 
formation. If the relation between optical surface brightness and gas density extends to still 
lower surface brightnesses this implies that galaxies that most easily escape detection in 
optical surveys are also the ones that might be missed in 21cm surveys. 

Until recently, blind H I surveys were not able to reach column densities much 
lower than lO^'^cm^^. These surveys would therefore still miss low density systems, and 
would not be capable to set strict constraints on the number density of LSB galaxies. 
However, the AHiSS is sensitive to column densities ~ lO^^cm"^ at a 5(7 level for gas 
filling the beam. Using an extrapolation of the scaling of Disney & Banks (1997), 
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A^Hi = 10^°(Mhi/-Lb)10°-^(^^~'*^\ where hb is the surface brightness averaged over the H I 
disk, and a typical value of Mm/ Lb — 1 for LSB galaxies, the AHiSS would be capable 
of finding galaxies with //^ = 32 magarcsec"^. Even if a ten times lower value of Mm/ Lb 
is used, and if the area of the galaxy were ten times smaller than the area covered by the 
Arecibo beam, still galaxies as faint as /^^ = 27 magarcsec"^ would be detectable. 

The VLA observations of the AHiSS galaxies can be used to set lower limits to average 
column densities for the sample. Due to the low spatial resolution of the observations 
('--' 1'). the measured column density is in many cases an average (A^hi) over the entire 
projected surface of the H I layer of the galaxy. If the galaxy is spatially unresolved, an 
upper limit to the extent Dm of the H I layer can be determined, leading to a lower limit 
of the average column density {Nm) oc Mhi-Djji • Although many of the galaxies in the 
sample are unresolved, and many of the true values of (A^hi) may be higher, we find no 
values of (Ahi) below 10^^-^cm~^. Hence, all galaxies in the sample have normal H I column 
densities, even though there are no observational selection effects against finding extreme 
low density systems. There is no indication of the existence of a group of extreme low 
column density galaxies that has been missed by previous H I surveys, simply because they 
were not capable of reaching the same low column densities as the AHlSS. Because we seem 
to be observing a lower limit to (A^hi), much higher than our detection limit, it is very 
unlikely that galaxies with even lower column densities are missed by the AHlSS. 

Theoretical predictions of the ionization of gas layers by the extragalactic UV 
background (Corbelli & Salpeter 1993, Maloney 1993, Charlton, Salpeter & hinder 1994) 
demonstrate physical mechanisms that can explain the possible non-existence of low column 
density neutral gas layers. These models predict photoionization by the extragalactic UV 
background of low column density regions, below lO^^'^cm"^. Furthermore, models by 
Quinn, Katz & Efstathiou (1996) show that the ionization only moderately suppresses the 
formation of galaxies with rotational speeds larger than 50 kms~^, but that it seriously 
affects the low density regions around these systems. The models arc verified by very deep 
VLA observations on one galaxy which appears to have a sharp truncation of the HI disk 
below a column density of 10^^'^cm^^ (van Gorkom et al. 1993). Further confirmation 
of ionization of low density HI comes from recent observations by Bland-Hawthorn et al. 
(1997) who have detected ionized gas beyond the truncated HI disk in NGC 253. 

5.6. HI self absorption 

The calculation of the total H I masses in this paper is based on the assumption 
that the optical depth of the H I layer is close to zero. Any possible influence of H I self 
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absorption, which will cause an underestimation of the true H I mass, is ignored. In this 
paragraph we will make a rough estimate of the influence of H I self absorption on the 
determination of the cosmological mass density f^Hi and the HiMF. The possible effect of 
self absorption could apply to all HiMFs compared in this paper (see Section |5.1| , ^ and 
Fig. 0), since none of these have addressed this possibility. 

The problem of self absorption for galaxies can be assessed statistically by plotting the 
21cm flux of different Hubble types as a function of inclination i to the line of sight. The 
line of sight through a inclined galaxy will be larger, generally causing a higher fraction 
of self absorption. Haynes & Giovanelli (1984) use data of a sample of 1500 galaxies with 
21cm fluxes measured with Arecibo. They define a correction factor /hi which is defined by 
the corrected flux divided by the measured flux and find a general trend: /hi = (cosi)"*^ 
where c is a constant dependent on morphological type. The values of c are found to be 
0.04 for Sa and Sab, 0.16 for Sb and 0.14 for Sbc and Sc galaxies. No correlations are found 
for morphological types earlier than Sa or later than Sc, indicating that self absorption 
is negligible in these types. Higher self absorption in types Sb to Sc can be explained by 
the fact that these galaxies generally have the highest H I surface densities (Cayatte et 
al. 1994). Furthermore, rotation curves of early type spiral galaxies show the strongest 
deviations from solid body rotation implying large velocity spreads per line of sight, which 
further decreases self absorption. 

Mean self absorption factors per morphological type can be obtained by averaging /hi 
over a random distribution of inclinations. This yields 



/o'^^^(cosi) sin idi 
Jq^'^ sin i di 



(/hi) = ■ . = V(l-c). (9) 



The mean correction factors then become 1.04 for Sa and Sab, 1.19 for Sb and 1.16 for Sbc 
and Sc galaxies. 

The cosmological H I mass density is dominated by high mass galaxies, which are 
statistically most likely to be late type spirals. The correction to Qm due to H I self 
absorption will therefore probably not be more than the value of (/hi) averaged over 
morphological types Sb to Sd. Assuming that all types Sb to Sd contribute equally to flm, 
we derive that the mean value of (/hi) is 1.10. Even in the most pessimistic approach, the 
correction factor can not be more than 1.19. 

What will be the effect of H I self absorption on the shape of the HiMF? Using the 
same arguments as above, we conclude that the effect on the high mass range will be 
marginal. The normalization factor 6* and the value which determines the kink, Mhi, will 
probably increase by no more than 10%. Galaxies that determine the faint end slope a are 
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mostly dwarf and LSB galaxies. On the one hand, self absorption may be unimportant 
because the gas density in these galaxies is usually low (van der Hulst 1993) and the dust 
content is presumably low (McGaugh 1994) which implies scarcity of clumped gas (Haynes 
& Giovanelli 1984). On the other hand, the rotation curves of dwarf and LSB galaxies 
often show solid body rotation (de Blok et al. 1996, Swaters 1997) which causes a low 
velocity spread along a line of sight, leading to high fractions of self absorption. These two 
counteracting effects make the value of (/hi) for low mass galaxies uncertain, but probably 
higher than that for high mass galaxies. After applying the correction factor to the low 
mass galaxies, some galaxies will shift to higher mass bins, eventually leading to a slightly 
shallower faint end slope of the mass function. The conclusion that the faint end slope of the 
mass function is shallow will therefore not be altered by the effects of H I self absorption. 

6. Conclusions 

We have used the Arecibo H I Strip Survey in combination with 21cm follow-up 
observations with the VLA to determine the H I mass function of gas rich galaxies in the 
local universe. The resulting HiMF is in good agreement with earlier estimates based on 
samples of optically selected galaxies. This implies that there is not a large population 
of gas rich dwarfs or low surface brightness galaxies, previously imdetected by optical 
surveys. The observed faint end slope of the HiMF is flat (a ~ 1.2) leaving no room 
for a large class of gas rich dwarfs. The cosmological mass density of H I in the local 
universe is ^lm{z = 0) = (2.0 ± 0.5) x also consistent with earlier estimates. The 

neutral gas content is dominated by high mass galaxies with lO^M© < Mm < lO^^M©. The 
observations indicate the existence of a lower limit to the average H I column density of 
19.7 cm~^, consistent with theoretical predictions concerning the ionizing extragalactic UV 
background. 

We thank G. Bothun, E. de Blok, P. Sackett, A. Szomoru, M. Verheijen, and the referee 
J. van Gorkom for useful comments. 
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TABLE 1 

Basic Parameters of Survey Galaxies. 



(1) 

Nr 


(2) 
Name 


(3) 
a (1950) 
[h m s) 


(4) 
* (1950) 

(° ' ") 


(5) 

Mai 
(log Mq) 


(6) 
Vq 
(km/s) 


(7) 
A* 

(") 


(8) 
Code 


Al 




19; 


:56: 


;40, 


.0 


14 


:07; 


:24 


s. 


.02 


1964 


33 




A2 




18: 


:37: 


;52. 


.6 


14 


:11; 


;53 


8. 


.32 


3991 


21 




A3 


UGC 11S20 


21; 


:47: 


;06, 


.0 


13 


:59; 


:50 


9, 


.12 


1107 


116 




A4-1 




21; 


:56: 


;41. 


.2 


13 


:59; 


:5S 


7, 


.99 


1720 


89 


A 


A4-2 


UGC 11S66 


21; 


:56: 


;09. 


.0 


13 


:53; 


:00 


9. 


.20 


1705 


507 


A 


A5 




20; 


:23: 


;45. 


.0 


14 


:06; 


:12 


9. 


.17 


4485 


49 




A6 


UGC 11611 


20; 


:3S: 


;11, 


.1 


14 


:05; 


:46 


9, 


.64 


5412 


47 




A7 


UGC 11652 


20; 


:55: 


;09, 


.1 


14 


:02; 


:5S 


9, 


.42 


4843 


52 




AS 


UGC 11617 


20; 


:41: 


;1S. 


.4 


14 


:07; 


:03 


9, 


.29 


5114 


152 




A9 


UGC 11921 


22; 


:06: 


;49. 


.2 


14 


:06; 


:53 


S. 


.82 


1672 


346 




AlO 


NGC 7437 


22; 


:55: 


;40, 


.9 


14 


:02; 


:2S 


S. 


.93 


2117 


151 




All 


UGC 1230S 


22; 


:5S: 


;49, 


.2 


14 


:04; 


:16 


9, 


.54 


2220 


261 




A12 




23; 


:23: 


;35, 


.4 


13 


:59; 


:09 


7, 


.92 


2860 


27 




A13 


UGC 11992 


22; 


:1S: 


;20. 


.7 


13 


:5S; 


:53 


9. 


.31 


3591 


115 




A14 


UGC 12705 


23; 


;34; 


;05. 


.2 


13 


:52; 


:48 


9. 


.67 


3968 


403 




A15 




00; 


:0S: 


;34. 


.1 


13 


:57; 


:45 


7, 


.35 


812 


100 




A16 




01; 


:30: 


;19, 


.2 


14 


:07; 


:05 


7, 


.68 


668 


392 




A17 




00; 


:17: 


;34. 


.1 


14 


:00; 


:50 


9, 


.18 


4787 


82 




AIS 




00; 


:21: 


;55. 


.2 


13 


:59; 


:05 


S. 


.51 


5397 


25 




A19 




00; 


:25: 


;2S. 


.S 


14 


:01; 


:30 


S. 


.90 


4552 


119 




A20 




00; 


:30: 


;51, 


.3 


13 


:59; 


:30 


S. 


.62 


5246 


3 




A21 




00; 


:41: 


;45. 


.9 


14 


:00; 


:51 


S. 


.33 


4372 


71 




A22 


UGC 1087 


01; 


:28: 


;46. 


.5 


14 


:01; 


;16 


9. 


.35 


4484 


45 




A23 




01; 


:39: 


;44. 


.4 


13 


:59; 


:20 


S. 


.79 


7103 


88 




A24 


NGC S20 


02; 


:05: 


;42. 


.3 


14 


:06; 


:46 


9, 


.56 


4418 


311 




A25 




02; 


:0S: 


;50, 


.7 


14 


:00; 


:17 


S. 


.62 


3794 


85 




A26 


UGC 1S17 


02; 


:1S: 


;47. 


.5 


13 


:5S; 


:22 


9. 


.63 


3735 


222 




A27 


UGC 2S39 


03; 


:41: 


;02, 


.0 


14 


:0S; 


:30 


9, 


.76 


6523 


156 




A2S 




04; 


:27: 


;13, 


.5 


14 


:09; 


:45 


S. 


.89 


4781 


68 




A29 




04; 


:30: 


;00, 


.0 


14 


:06; 


:40 


9, 


.00 


4523 


118 




A30 


UGC 1294 


01; 


:47: 


;26, 


.2 


22 


:54; 


:44 


S. 


.71 


2861 


73 




A31 




01; 


;50; 


;57. 


.5 


22 


:53; 


;59 


8. 


.97 


4956 


124 




A32 


F 477-01 


01; 


:51: 


;46. 


.9 


22 


:57; 


:35 


9, 


.13 


4989 


90 




A33 


UGC 193S 


02; 


:25: 


;32, 


.1 


22 


:59; 


:20 


9, 


.48 


6383 


108 




A34 


UGC 2927 


03; 


:5S: 


;41. 


.5 


22 


:5S; 


:20 


9, 


.81 


6251 


180 




A35 




04; 


:11: 


;36. 


.0 


23 


:00; 


:00 


9. 


.21 


5464 


148 




A36 




04; 


:44: 


;14. 


.3 


23 


:04; 


:19 


8. 


.45 


5088 


1 




A37-1 




04; 


:47: 


;54. 


.5 


23 


:10; 


:50 


S. 


.98 


4272 


375 


B 


A37-2 


UGC 31S3 


04; 


:4S: 


;1S. 


.6 


23 


:07; 


:29 


S. 


.93 


4390 


192 


B 


A3S 




04; 


:49: 


;52. 


.4 


23 


:05; 


:12 


9. 


.22 


4416 


33 




A39 




04; 


:59: 


;39, 


.3 


23 


:10; 


:10 


9, 


.43 


6339 


294 




A40 




05; 


:52: 


;39, 


.4 


23 


:02; 


:29 


9, 


.06 


5812 


367 




A41 




06; 


:11: 


;31, 


.1 


23 


:04; 


:49 


9, 


.63 


5902 


299 




A42 




06; 


:3S: 


;5S, 


.9 


23 


:12; 


:05 


S. 


.38 


7012 


32 




A43 


UGC 3751 


07; 


;10; 


;53. 


.4 


23 


:10; 


;05 


8. 


.94 


2300 


208 




A44 


CGCG116-001 


06; 


:50: 


;06, 


.2 


23 


:11; 


:30 


S. 


.11 


4575 


40 




A45-1 




06; 


;51; 


;09, 


.6 


23 


:11; 


;52 


S. 


.32 


4593 


28 


C 


A45-2 




06; 


:51: 


;26, 


.4 


23 


:10; 


:3S 


s. 


.63 


4562 


101 


C 


A46 




07; 


:07: 


;27. 


.5 


23 


:13; 


:02 


s. 


.90 


5457 


18 




A47 




07; 


:11: 


;33, 


.6 


23 


:14; 


:49 


9, 


.00 


5473 


74 




A48 




07; 


:43: 


;23, 


.7 


23 


:12; 


:56 


9, 


.22 


5491 


151 




A49 


UGC 4031 


07; 


:45: 


;19, 


.9 


23 


:21; 


:51 


9, 


.29 


6834 


378 




A50 




OS; 


:06: 


;52, 


.0 


23 


:17; 


:47 


S. 


.41 


4790 


61 




A51 




08; 


;10; 


;56. 


.9 


23 


:17; 


:44 


8. 


.01 


4667 


60 




A52 


UGC 4299 


OS; 


:13: 


;00, 


.2 


23 


:21; 


:37 


9, 


.31 


4107 


271 


D 


A53 




OS; 


:13: 


;29, 


.0 


23 


:17; 


:51 


S. 


.61 


4143 


44 


D 


A54 




OS; 


:20: 


;1S. 


.6 


23 


:1S; 


:29 


S. 


.29 


4029 


60 




A55 


UGC 4405 


08; 


;23; 


;36. 


.2 


23 


:21; 


;25 


9. 


.49 


5617 


276 




A56 


CGCG120-049 


OS; 


:51: 


;26, 


.S 


23 


:20; 


:55 


S. 


.15 


2882 


114 




A57 




10; 


:31: 


;55, 


.5 


23 


:20; 


:54 


S. 


.00 


1238 


101 




ASS 




09; 


:33: 


;37, 


.4 


23 


:20; 


:03 


S. 


.81 


6132 


45 




A59-1 


NGC 2929 


09; 


:34: 


;39. 


.3 


23 


:23; 


:1S 


9. 


.72 


7499 


148 


E 


A59-2 


NGC 2931 


09; 


:34: 


;46. 


.S 


23 


:2S; 


:06 


9, 


.45 


7478 


436 


E 


A59-3 


NGC 2930 


09; 


:34: 


;40, 


.S 


23 


:26; 


:12 


9, 


.26 


7376 


321 


E 


A60 


UGC 5214 


09; 


:42: 


;32, 


.5 


23 


:1S; 


:07 


9, 


.39 


5480 


180 




A61 


UGC 549S 


10; 


:09: 


;16 


.4 


23 


:19; 


:42 


9. 


.14 


6306 


138 
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